Emergence of Secondary Motifs in Tube-Like Polymers in a Solvent 
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We study the effects of two kinds of interactions in tube-like polymers and demonstrate that they 
result in the formation of secondary motifs. The first has an entropic origin and is a measure of the 
eflFective space available to the solvent. The second arises from solvophobic interactions of the solvent 
with the polymers and leads to an energy proportional to the contact surface between the tube and 
solvent particles. The solvent molecules are modeled as hard spheres and the two interactions 
are considered separately with the solvent density affecting their relative strength. In addition to 
analytical calculations, we present the results of numerical simulations in order to understand the 
role played by the finite length of short polymers and the discrete versus continuum descriptions of 
the system in determining the preferred conformation. 

PACS numbers: 87.14.Ee, 05.70.-a, 61.20.-p 
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The study of entropic effects has a venerable history in 
physics. Onsager JL] studied the nature of an entropically 
driven isotropic-nematic transition of a fluid composed of 
infinitely thin hard rods. The addition of a small amount 
of polymers to a colloidal suspension induces an effective 
attraction between the colloidal particles which can lead 
to flocculation through the Asakura-Oosawa (AO) mech- 
anism Our focus, in this letter, is on studying the 
optimal conformations of a tube subject to compaction 
through its interactions with the solvent. Our principal 
results are that the tube adopts helical and planar con- 
formations with the latter resulting from the discreteness 
of the tube. 

Our results are useful for understanding proteins [1, 01 ■ 
In spite of their complexity, proteins fold into a lim- 
ited number Q of evolutionarily conserved structures 
made up of helices and almost planar sheets. It has 
been proposed that the geometries of protein folds orig- 
inate from the common attributes of all proteins that 
can be encoded in a simple geometrical model of a flex- 
ible tube of non-zero thickness 0, B 0- Proteins oc- 
cupy the marginally compact phase of a tube-like poly- 
mer [l^, [ll|, [12, [l^ which is rich in secondary mo- 
tifs and has a relatively low structural degeneracy com- 
pared to the generic compact polymer phase. While it is 
well recognized that the formation of secondary structure 
arises from backbone-backbone hydrogen bonds [3, [l^ , 
the solvent is known to play a fundamental role in pro- 
moting the folding process The importance of the 
solvation free energy and of entropy induced interactions 
has been increasingly recognized in the last few years 

In a system of hard spheres, the excluded volume con- 
straint can be enforced by forbidding configurations in 
which the distance between pairs of sphere centers is 
smaller than the sphere diameter. In the continuum 



limit, a self interacting tube- like polymer cannot have 
two body interactions ^] . The self-avoidance of a flexi- 
ble tube of thickness A (the tube radiusjcan be enforced 
through a three-body potential 0, HO, Hlj which forbids 
conformations in which the radius of a circle through 
an arbitrary triplet of points on the axis is smaller than 
the tube thickness. A solvophobic polymer in a solvent, 
such as a protein in vivo, adopts a self-avoiding confor- 
mation in which the contact surface is minimized, an ef- 
fect that may be captured by considering Q an effective 
solvent induced (free-)energy proportional to the contact 
surface between the tube and the solvent particles. Fur- 
thermore, due to the discrete nature of the solvent one 
expects ent ropy -induced effective interactions of the AO 
type 0, [H, [l3l as investigated by Snir and Kamien [l3| 
and more recently by Hansen-Goos et al. 19] for a tube- 
like polymer in helical conformations. The AO type of 
interaction yields an effective (free-) energy proportional 
to the volume excluded to the solvent. 

It was noted recently by Hansen-Goos et al. [l^, [l^l 
that Hadwiger's theorem of integral geometry [2^, under 
rather general conditions, allows one to write the solva- 
tion free energy Fsoi of the tube in a solvent as 



Fsoi = PV,, 



kC + kX 



(1) 



where Vexc, ^acc, C and X are the volume excluded to 
the solvent, the area accessible to the solvent and inte- 
grated mean and Gaussian curvatures of the accessible 
surface, respectively. P and a are the solvent pressure 
and the planar surface tension respectively, whereas k 
and k represents the two bending rigidities; the proper- 
ties of the solvent enters through these coefficients and 
can hence be computed separately. The AO approxima- 
tion amounts to setting a — k — k = and assuming an 
ideal gas solvent P = Pid corresponding to an infinitely 
diluted system. The principal aim of our letter is to pro- 



2 



vide exact expressions for the excluded volume and the 
contact surface of a tube- like polymer in a solvent of hard 
spheres and to investigate the effects of these terms both 
analytically and by Monte Carlo simulations in determin- 
ing the optimal conformations of short polymers. 

In a continuum description, the axis of a tube is de- 
scribed by the curve R(s) parameterized by its arc length 
s. The Frenet reference frame at location s is given in 
terms of the unit vectors {T(s), N(s), B(s)}, the tan- 
gent, the normal and the binormal respectively The 
curvature at location s, k{s) (the inverse radius of curva- 
ture), is defined in terms of one of the Frenet-Serret equa- 
tions, dT (s) /ds = K (s) N (s) and B (s) = T (s) x N (s). 

The position vector of a point inside the tube of length 
L is 



R(s) + p cos6' N(s) -t-sin6l B(s) (2) 



where < 6* < 27r, < p < A and < s < L. 
The volume element in the above curvilinear coordi- 
nate system is given by dV = p\x (s, p, 0) \ dsdpdO, where 
X (s, p,9) = 1 — k{s)p cosd. Let e be the radius of the sol- 
vent particles assumed to be spherical. For a straight 
tube the excluded volume is simply Vexc,s = t^LA^, 
where = A -|- e, since the volume inaccessible to the 
centers of the solvent particles is a tube with thickness 
Aj. In order to calculate the excluded volume when the 
tube curls and different regions of the tube come close to 
each other, we have to consider points which are shared 
by more than one circular section of radius A^ associated 
with distinct points of the tube axis. We avoid multi- 
ple counting by introducing the integer function n{s, p, 9) 
which gives the number of circular cross sections of radius 
Aj, centered at distinct points of the tube axis, contain- 
ing the point r{s,p,6) where p < Ag. A point, r{s,p,6), 
belongs to the circular section associated with R(s') if 
the vector distance ds,p,e{s') = r{s,p,9) — R(s') is per- 
pendicular to the tangent at s', T(s'), and its modulus is 
less than A^. Note that n{s, p,9) > 1 because s' = s triv- 
ially satisfies the above conditions. The excluded volume 
is thus given by 



V n{s,p,0) 



(3) 



where T) is the domain < < 27r, < /? < A^ and 
< s < L. The same integer function can be used to 
obtain an expression for the contact surface on setting 
the p coordinate equal to the radius A of the tube. 



Scont = / dV 6{p 
I'D 



A)e[l-n(s,A„0)] (4) 



Here Q{x) is the step function equal to 1 if a; > and 
otherwise. It is not difficult to show that Eq. (g]) is 
equivalent to the formula presented in Ref . , for both 
an infinite length tube and a finite one. In the latter 
case hemispheres needs to be attached to both ends while 
computing n{s,p,9) > 1, as is the case in our numerical 



simulations. In a swollen configuration corresponding to 
the absence of overlap, n{s,p,9) = 1 for all {s,p,9) and 
k(s) < A7^ for all s, and Eqs. ^ and ^ reproduce the 
correct results Vexc,s = ttXA^ and 'Scant, s — 27rAL. It is 
illuminating to rewrite the difference Vexc.s ~ Vexc as 

AV^ = Vexc.s - Vexc = Voy " 2K,<0 (5) 

where we have introduced the overlap volume 

n{s,p,9) - 1 



V n{s,p,9) 



(6) 



which is the part of the excluded volume with non- 
vanishing overlap, i.e. with n{s,p,9) > 1, and the vol- 
ume in the highly curved region, i.e. k{s) > A^T^, 



dv e[-xis,p,9)] 



(7) 



As a result, the difference AV between the excluded vol- 
umes in a swollen and a compact conformation is not 
simply equal to the overlap volume Vov as in the usual 
Asakura-Oosawa mechanism involving rigid spheres j^]- 
There is no change in shape of the colloidal particle 
whereas the tube polymer undergoes a shape-change 
upon folding. This difference has some interesting conse- 
quences. Suppose, for instance, that the tube is twisted 
and connected at the two ends to form a torus/donut of 
thickness A and radius R. If i? > A^ one clearly has 
Ay = 0. Ifi?<Ae, a straightforward calculation shows 
that Vov = V^<o so that AV" = -$(i?/Ae)47rA3/3 be- 
comes negative with the function 



1 



-X arccosx 



(8) 



positive in the interval < a; < 1. The swollen con- 
figuration is entropically more favorable compared to a 
tight donut. This counterintuitive result stems from the 
fact that the curvature yields both inflated and deflated 
volume elements. In the absence of overlap, the two ef- 
fects balance against each other yielding AV ~ 0. When 
X < 0, the same balance can be achieved only by allow- 
ing deflated "negative" volumes. An imbalanced inflated 
volume is then present, which may overturn the excluded 
volume reduction due to the presence of overlap, as in the 
donut case. 

Consider an infinite helix of radius R and pitch P. 
The translational invariance along the tube axis simpli- 
fies the calculations, and we can consider both the ex- 
cluded volume and the contact surface per unit length. 
Fig. [T] shows the values of c = P/R for the helices 
which minimizes Vexc/L^ as the solvent molecule ra- 
dius changes - the analogous results for Y,cont/ L exhibits 
similar behavior. When e is below a threshold value, 
the optimal co nfig uration corresponds to the ideal helix 
0j Elj IitI . [TqI . |25| | which simultaneously minimizes the 
local and non-local radius of curvature 2l| . For an ideal 
helix, c = c* = 2.5126 .. . which is within 3% of the value 
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FIG. 1: Plot of the pitch-to-radius value c — P/R for the op- 
timal helix on minimizing the excluded volume as a function 
of the ratio e/A, for different helix lengths, L. Below a given 
threshold value (shown for the L — oo case) of the solvent 
radius, the ideal helix conformation is obtained. Our results 
were obtained by first setting the value of the pitch, P, so 
that two different turns of the helix have minimum distance 
for a fixed 7? Second, we identify the value of c which 

minimizes Vexc/L for a given e and these values are plot- 
ted against each other. An analogous calculation for T.cont/L 
shows a similar trend towards optimal finite helices having a 
greater c with respect to the infinite helix. In the inset, the 
difi'erent behaviors obtained on considering the minimization 
of the excluded volume and the maximization of the overlap 
volume is shown. 



for a— helices in proteins Q. Interestingly, a similar op- 
timality condition is also found for double strand DNA 
[261] . These observations suggest that a— helices might 
originate in general optimization processes transcending 
the large differences in the amino acid sequences of dif- 
ferent proteins. As the solvent sphere size e increases 
beyond the threshold. Fig. [T] shows that the optimal 
value of c decreases from its value in the ideal helical 
structure. Our results confirm those obtained in previous 
work [13, [1^ . For the the excluded volume interaction, 
the helix maintains its ideal pitch-to-radius ratio c* up 
to a critical value e* ~ 0.0835A. Above this threshold, 
the ideal helix is no longer the optimal conformation - 
even though the excluded volume along the central axis 
increases, there is a better overlap between subsequent 
turns thus decreasing the excluded volume. The inset of 
Fig. [T] shows the results maximizing the overlap rather 
than minimizing the excluded volume. The optimal value 
c* is maintained up to e* ~ O.IA, a value larger than that 
obtained in the minimum excluded volume case. This 
mirrors the situation with the donut: for large curva- 
ture, there is a non-zero, albeit small, difference between 
AV and Vov as dictated by Eq. An exact analytical 
computation for the contact surface [2^ yields an iden- 
tical pattern with a smaller critical value e* = 0.04627A 



for the solvent radius. Again, this is in perfect agreement 
with the value given in Ref. ^19.] using a different method. 
Our method allows one to treat the case of finite helices 
efficiently using numerical integration. Fig. [1] shows the 
optimal value of c for helices of varying length, L. Short 
helices show a tendency to remain ideal for higher values 
of the solvent radius. 

When the tube is very long or there are many inde- 
pendent tubes, one might expect that the helical confor- 
mation is destabilized. For many long tubes, an optimal 
configuration which minimizes both the contact surface 
and the excluded volume is likely a regular arrangement 
of parallel tubes forming a hexagonal arrangement some- 
what reminiscent of the Abrikosov phase in Type II su- 
perconductors [iH. We have verified that > 4 in- 
finitely long tubes have a lower excluded volume per unit 
length with respect to the optimal helix (for a single tube 
one should also take into account the — 1 turns) . 
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FIG. 2: Conformations adopted by tubes of length (a) L = 
lOA and (b) L = 20A subject to either the excluded volume or 
contact surface prescriptions for promoting compaction. The 
values of e/A are 0.005, 0.025, 0.05, 0.10, 0.20 and 0.50. The 
resulting conformations for the shorter tube are either saddles 
or helices for e/A > 0.20, while for smaller e, after a crossover 
phase, the hairpin becomes the ground state - the planarity of 
the structure is highlighted on the left of (a). For the longer 
tube, the planar phase consist of /3-sheet-like structures and 
the helicoidal phase is characterized by double helices, saddles 
and irregular helices. The planar phase disappears for the 
contact surface case. All these conformations are akin to those 
found in Ref. [ll. 



4 





L = 


lOA 


L = 


20A 


e/A 


AV/(LA2) 


AE/(LA) 


AV/(LA2) 


AE/(LA) 


0.005 


0.049 


0.62 


0.053 


0.98 


0.025 


0.061 


0.80 


0.067 


1.2 


0.050 


0.079 


0.98 


0.098 


1.5 


0.100 


0.12 


1.2 


0.18 


1.9 


0.200 


0.24 


1.6 


0.39 


2.1 


0.500 


0.77 


2.1 


1.3 


2.9 



TABLE I: Values of AV"/(LA2) and AE/(LA) related to the 
optimal structures, for the tested values of e/A and L/A. 

In order to investigate the possibility for other kinds 
of optimal conformations for a short tube as well as the 
role of discrete versus continuum descriptions of the tube, 
we resort to Monte Carlo simulations using the formu- 
lae given above for both the excluded volume and the 
contact surface. We seek to identify the conformations 
which minimize Vexc or Y^cont as the solvent radius e and 
the polymer length L are varied (the thickness A is kept 
fixed because the results depend only on the dimension- 
less ratios e/A and L/A). Our Monte Carlo simulations 
combine crankshaft and pivot moves of a homopolymer 
with a Metropolis annealing schedule [1^ following ini- 
tializing the chain in an arbitrary extended conforma- 
tion. We have considered two different polymer lengths, 
L = lOA and L = 20A, fixing the bond length b = A/2. 
The results are summarized in Fig. [2]- the optimal struc- 
tures obtained on varying e are shown ~ while the values 
of AV^ = Vexc.s ~ Vexc and AE = T.cont,s - Scont reached 
are reported in Tab. [D Fig. [5] (a) displays the conforma- 
tions obtained with L — 10 A. For both the excluded vol- 
ume and contact surface, two different "phases" emerge: 



at small e, a planar phase, in which the ground state 
has a typical planar hairpin structure, is found; at large 
e, a helicoidal phase is found in which helices and sad- 
dles dominate, with comparable Vexc and 'Scant- The 
two different regimes are separated by a crossover re- 
gion characterized by the coexistence of all these struc- 
tures. Note that the excluded volume effect exhibits a 
stronger propensity for the planar phase. This is evident 
for L — 20A shown in Fig[2](b) - the contact surface pla- 
nar phase disappears already at e/A = 0.005, the lowest 
value tested. The helicoidal phase in this case consist of 
double helices, saddles and irregular helices, with turns of 
different lengths. Our simulations suggest the existence 
of an energy barrier between the two classes of confor- 
mations (hairpin versus helix or saddle) in the crossover 
region. The appearance of planar structures is a con- 
sequence of the discrete nature of the polymer (non-zero 
bond length) which plays a crucial role especially at small 
e. 

In summary, tube-like polymers with naturally arising 
solvent induced interactions exhibit low free energy con- 
formations with secondary motifs. This suggests that the 
secondary motifs commonly found in bio-polymers such 
as proteins and DNA have a common and fundamental 
origin which transcends chemical detail. When the dis- 
crete nature of the polymer dominates, planar conforma- 
tions akin to the /9-sheet in proteins emerge along with 
helical conformations as the optimal ones. In contrast, 
when the solvent molecule size is sufficiently large so that 
the discrete nature of the polymer can be neglected, the 
continuum approximation is valid and helical conforma- 
tions dominate with single and double helix conforma- 
tions. 
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